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Abstract. The recently introduced GALI method is used for rapidly detecting 
chaos, determining the dimensionality of regular motion and predicting slow dif- 
fusion in multi-dimensional Hamiltonian systems. We propose an efficient com- 
putation of the GALIfe indices, which represent volume elements of k randomly 
chosen deviation vectors from a given orbit, based on the Singular Value Decom- 
position (SVD) algorithm. We obtain theoretically and verify numerically asymp- 
totic estimates of GALIs long-time behavior in the case of regular orbits lying 
on low-dimensional tori. The GALI^ indices are applied to rapidly detect chaotic 
oscillations, identify low-dimensional tori of Fermi-Pasta-Ulam (FPU) lattices 
at low energies and predict weak diffusion away from quasiperiodic motion, long 
before it is actually observed in the oscillations. 

1 Introduction 

A great variety of physical systems can be described by Hamiltonian systems or symplectic 
maps |H2l3j . Their applications range from the stability of the solar system [4] and the con- 
tainment of charged particles in high intensity magnetic fields [3] to the blow-up of hadron 
beams in high energy accelerators [5] and the study of simple molecules and hydrogen-bonded 
systems |6I7) . More recently, much attention has been focused on the dynamics of nonlinear 
lattices and in particular on localization in the form of q-breathers energy transport and 
cquipartition properties of Fermi-Pasta-Ulam (FPU) particle chains [9]. The main difficulty 
with determining the nature of the motion in Hamiltonian dynamics is that regular and chaotic 
orbits are distributed in phase space in very intricate ways, in contrast with dissipative systems, 
where all orbits eventually fall on regular or chaotic attractors. 

The most widely used method for distinguishing order from chaos in dynamical systems is 
the evaluation of the Lyapunov Exponents (LEs) <Ji for each given orbit. Benettin et al. [TO] 
studied the problem of determining all LEs theoretically and proposed in [llj an algorithm for 
their numerical computation. In practice, one often calculates only the largest LE, a\, following 
one deviation vector from the given orbit and if a% > the orbit is characterized as chaotic. 

In the present paper, we apply the Generalized Alignment Index of order k, GALIfc, k = 
2, . . . , 2N [T2], to distinguish efficiently between regular and chaotic orbits of multi-dimensional 
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Hamiltonian systems of N degrees of freedom. These indices generalize a similar indicator 
called Smaller ALignment Index (SALI) |13j . in that they use information of more than two 
deviation vectors from the reference orbit (see [12] for more details). In particular, the GALIfc 
is proportional to volume elements formed by k initially linearly independent unit deviation 
vectors whose magnitude is normalized to unity at every time step. 

After recalling briefly the definition of GALIfc and its general behavior for regular and 
chaotic motion in Sect. 12. 1[ we propose in Sect. 12.21 a technique for the efficient computation 
of GALIfc, based on the Singular Value Decomposition (SVD) of the matrix having as rows 
the k normalized deviation vectors. In Sect. 12.31 we study theoretically the behavior of GALIs 
for regular orbits of N degrees of freedom Hamiltonian systems that lie on an s-dimensional 
torus, with 1 < s < N. In Sect. [3J we present applications of the GALIfc approach to an 
N = 8 particles FPU lattice and show that the GALIfc indices: (i) can be used for the rapid 
discrimination between regular and chaotic motion, (ii) determine the correct dimensionality 
of tori and (iii) predict that the motion is weakly chaotic, long before this can be seen in the 
actual oscillations. 

By 'weak chaos' we mean weakly diffusive motion through a network of resonances. On the 
other hand, when one speaks of 'sticky' orbits [14j . one generally refers to motion occurring 
just outside the boundary of a large regular region, where orbits remain for a long times before 
rapidly escaping to distant parts of phase space through a large chaotic sea. In our case, 
'weakly chaotic' motion occurs within regimes 'surrounded' by (often high-dimensional) tori, 
is characterized by very small LEs and is reminiscent of what is called in the literature Arnol'd 
diffusion [3J. Finally, in Sect.[3J we summarize the results and present our conclusions. 

2 The GALI method 

2.1 Definition and behavior of GALI 

Following [12] let us first briefly recall the definition of GALI and its behavior for regular and 
chaotic motion. We consider a Hamiltonian system of N degrees of freedom having a Hamilto- 
nian H(qi, . . . , (iNiPit ■ ■ ■ tPn) where qi and pi, i = 1, . . . , N are the generalized coordinates and 
momenta respectively. An orbit of this system is defined by a vector x(t) = (xi(t), . . . , X2N{t)), 
with Xi = qi, Xi+N = Pi, i = 1, ...,iV. This orbit is a solution of Hamilton's equations 
dx/dt — V(x) = (dH/dp , — dH/dq), while the evolution of a deviation vector w(i) from 
x(i) obeys the variational equations dw/dt = M(x(i))w, where M = dV/dx is the Jacobian 
matrix of V. 

Let us follow k normalized deviation vectors wi, . . ., Wk (with 2 < k < 2N) in time, and 
determine whether they become linearly dependent, by checking if the volume of the correspond- 
ing fc-parallelogram goes to zero. This volume is equal to the norm of the wedge or exterior 
product [15j of these vectors. Hence we are led to define the following 'volume' element: 



which we call the Generalized Alignment Index (GALI) of order k. We note that the hat ( A ) 
over a vector denotes that it is of unit magnitude. Thus, for each initial condition x(0), we solve 
Hamilton's equations for x(t), together with their variational equations for k initially linearly 
independent deviation vectors Wi, i = 1, . . . , k. Clearly, if at least two of these vectors become 
linearly dependent, the wedge product in ([1]) becomes zero and the volume element vanishes. 

In the case of a chaotic orbit all deviation vectors tend to become linearly dependent, aligning 
in the direction of the eigenvector which corresponds to the maximal Lyapunov exponent and 
GALIfc tends to zero exponentially following the law [12"] : 



GALIfc(t) = ||i&i(t) Aii 2 (i) A--.At& fc (t)|| , 



(1) 



GALI fc (t) 



e 



[(fl— <T2) + (f 1—<T3) H h(fl <7fc)]t 



(2) 



where a~x , . . . , o~k are approximations of the first k largest Lyapunov exponents. In the case 
of regular motion on the other hand, all deviation vectors tend to fall on the TV-dimensional 
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tangent space of the torus on which the motion lies. Thus, if we start with k < N general 
deviation vectors they will remain linearly independent on the A-dimensional tangent space 
of the torus, since there is no particular reason for them to become aligned. As a consequence 
GALIfe remains practically constant for k < N. On the other hand, GALIfe tends to zero for 
k > N, since some deviation vectors will eventually become linearly dependent, following a 
particular power law which depends on the dimensionality N of the torus and the number k of 
deviation vectors. So, the generic behavior of GALIfe for regular orbits lying on ^-dimensional 
tori is given by |12j : 

f constant if 2 < k < N 
GALIfe (t) ~ | . f N < k < m . (3) 



2.2 Numerical computation of GALI 

In order to numerically compute GALIfe in [TJ] using ([T]) we considered as a basis of the 
2Af -dimensional tangent space of the Hamiltonian flow the usual set of orthonormal vectors 
ei = (1, 0, 0, . . . , 0), e 2 = (0, 1, 0, . . . , 0), &2N = (0, 0, 0, . . . , 1). So, any unitary deviation 
vector Wi can be written as: 



2N 
.7=1 



» = 1,2,. 



(4) 



where Wij are real numbers satisfying Y^ij=i w fj = 1- Considering the k x 2A^ matrix A having 
as rows the coordinates of k such vectors, we can write equations Q in matrix form as: 



~Wx 













wn W12 ■■■ W 12 N 

W21 W22 ' ■ ■ W22N 
Wkl W k2 ■■■ W k2 N 





ei 




ei 




e-2 




e-2 






= A 






-&2N - 







(5) 



The norm of the wedge product of the k deviation vectors appearing in {T]) was defined in [12] 
as the square root of the sum of the squares of the determinants of all possible k x k submatrices 
of matrix A. So, for the computation of GALIfe we have: 



GALIfe = llt&i A t&2 A • • • A «)fc|| = { ^2 

l<ii<i 2 <-<ik<2N 



det 



V 



win w ii 2 
W2i x w 2i 2 

Wkn w ki2 



w 2lk 



1/2 



. (6) 



where the sum is performed over all possible combinations of k indices out of 2N . 

Eq. ([6]) is ideal for the theoretical determination of the asymptotic (long time) behavior 
of GALIs for chaotic and regular orbits (see [T^] for more details), as well as, for regular 
orbits that lie on low dimensional tori (see Sect. 12.31 below). However, from a practical point 
of view the computation of determinants is not the most efficient way of computing GALIfe, 
as has already pointed out in [16j . For low dimensional systems, the number of determinants 
appearing in ((6|) for the computation of GALIfe is not prohibitive, but as the number of degrees 
of freedom increases the computation can become extremely time consuming and in some cases 
impractical. For example, in the case of an A^ = 8 degree of freedom system, like the one studied 
in Sect. [3l the computation of GALIg requires the evaluation of 12870 8x8 determinants, while, 
GALI15 in an A^ = 15 degree of freedom system requires the computation of 155117520 15 x 15 
determinants! 

We have already mentioned that GALIfe measures the volume of the fc-parallelogram P k 
having as edges the k unitary deviation vectors Wi, i — l,...,fc of ((4]) and ([5|) above. The 
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volume of Pk is then given by (see for instance [IT]): 



vol(P fe ) = ^/det(A-A 1 ) , (7) 

where ( T ) denotes transpose. Since det(A • A T ) is equal to the sum appearing in (J6j) (this 
equality is called Lagrange's identity, see for instance [18]), we have: 



GALIfc = y det(A • A T ), (8) 

as an alternative way of computing GALIfc , where only the multiplication of two matrices and 
the square root of one determinant appears. 

A different way of evaluating GALIfc, which actually proved to be more accurate, is obtained 
by performing the Singular Value Decomposition (SVD) of A T . So, the 2N x k matrix A T can 
be written as the product of a 2N x k column-orthogonal matrix U, a k x k diagonal matrix Z 
with positive or zero elements Zi, i = 1, . . . , k (the so-called singular values), and the transpose 
of a k x k orthogonal matrix V: 

A T =U-Z-V T . (9) 

We note that matrices U and V are orthogonal so that: 

U T -U = V T -V = Ifc, (10) 

with Ifc being the k x k unit matrix. For a more detailed description of the SVD method, as 
well as an algorithm for its implementation the reader is referred to [19] and references therein. 
Using Eq. ([5]) for the computation of GALIfc, as well as Eqs. © and (fit)]) , we get: 



GALL. = \ / det (A • A T ) = J det ( V- Z T • U T U Z V 



det (V • diag(^) • V T ) = v / det(diag(z l 2 )) = f[ 



Thus, we conclude that GALIfc is equal to the product of the singular values of matrix A ([5]) 
defined by the k normalized deviation vectors and at the same time, we theoretically explain 
the computationally verified equality of this product to values of GALIfc computed by © which 
was reported in [16] . The SVD approach, therefore, provides a very accurate determination of 
the logarithm of GALIfc, which can now be used for the discrimination between regular and 
chaotic motion, as it leads to: 

k 

log(GALIfc) = 5>g(z,) , (11) 
»=i 

and for this reason we implement this approach for the computation of GALIfc in numerical 
applications. From Eq. ([lip we see that the computation of the singular values Zi with the 
usual double precision accuracy, permits the accurate determination of large negative values of 
log(GALIfc), which correspond to very small values of GALIfc. 

The problem of the numerical computation of orbits, keeping constant the numerical value 
of the Hamiltonian H, is of great importance in numerical studies of dynamical systems. Several 
integration schemes have been developed and applied over the years with varying degrees of 
success (see [20] for a survey of such methods). In our study we use an 8th order Runge- 
Kutta method proposed in [21] , both for the numerical integration of Hamilton's equations 
(evolution of an orbit), as well as for the variational equations (evolution of deviation vectors). 
The scheme proved to be very efficient and accurate since, in all our computations, we always 
kept the relative error of the values of the Hamiltonian function (\H(i) — H(0)\/\H(0)\) below 

io- 12 . 
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2.3 Behavior of GALI for regular orbits of low dimensional tori 

Now we turn to the GALI& method for regular orbits of an N degree of freedom Hamiltonian, 
that lie on s-dimensional tori, with 1 < s < N. In this case, one could perform a local transfor- 
mation to action-angle variables, Ji,#i, whence Hamilton's equations of motion can be easily 
integrated to give Ji(t) = Jjo, 0i(t) = 6io + u>i( Jio, . . . , J s o) t, where Jjo, Oio are the initial 
conditions (i = 1, . . . , N) and u>i = for s < i < N. Denoting by £j, r\i small deviations from 
the Ji and 9i respectively and using as basis of the 2A-dimensional tangent space of the Hamil- 
tonian flow the 2N unit vectors {i>i, . . . , V2n}, such that the first TV of them correspond to the 
N action variables and the remaining ones to the N conjugate angle variables, any deviation 
vector w,-, i = 1, 2, . . . can be written as 

N N / N \ 

= +E ? ?K°) + 5>*&(°)* (12) 

3 = 1 J = l \ k=l I 

with LUkj = duJk/dJj computed from the initial values Jjo for k,j = 1,2, ... ,s and uikj = 
for k,j = s + l,s + 2, . . . ,N. It follows easily from the above that for sufficiently long times, 
||w(f)||~t. 

Let us now study the case of fc, initially linearly independent, randomly chosen, unit 
deviation vectors {t&i, . . . , Wk} expressed in terms of the new basis by the transformation 
[u>i . . . Wk ] T = D ■ [ii . . . V2n] T - The random choice of the initial deviation vectors corre- 
sponds to the generic (most probable) configuration that none of them lies in the tangent space 
of the torus. In the opposite case, the results do not change qualitatively since GALL still 
exhibit power law decays, but with slightly different exponents [T2]. Defining by £ k and r\\ the 
k X 1 column matrices of initial conditions, the matrix D assumes the form 

D(t) ~ i • T>\t) = i [i\ . . . i% Vl + E"=i ■■■V L s+ ELi V h s+ i ■ ■ ■ V k N ] - 

(13) 

where we have replaced the first factor on the right by its asymptotic expression for long times. 

In the case of an s-dimensional torus, the k deviation vectors eventually fall on its s— 
dimensional tangent space spanned by vn+i, ■ ■ ■ ,vn+s- If we start with 2 < k < s deviation 
vectors, since there is no reason for them to become linearly dependent, their wedge product 
yields GALIfe indices that are different from zero. However, if we start with s < k < 2N 
deviation vectors, some of them will necessarily become linearly dependent and thus their 
wedge product (as well as the GALL-) will tend to zero not exponentially but by a power law, 
as we explain below. 

In order to determine the leading order behavior of the GALL-, we search for the fastest 
increasing determinants of all k x k minors of the matrix D fe , as t grows. For 2 < k < s, these 
determinants have k columns chosen among the s columns of matrix T) k with u>ij ^ and grow 
as t k , thus providing constant terms to the GALIfc. All other determinants contain at least one 
column from the 2N — s time independent columns of matrix D fc and introduce terms that 
grow slower than t k , having ultimately no bearing on the behavior of GALIfe(t). This yields the 
important result that GALIfc(i) ~ constant for 2 < k < s. 

Next, we turn to the case of s < k < 2N — s. The fastest growing determinants are again 
those containing the s columns of the matrix D fe with wy 0. The remaining k — s columns 
are chosen among the 2(N — s) columns of D fe which are time independent (excluding the £^ 
columns with i < s). Among these determinants, the fastest increasing ones are those containing 
as many columns proportional to t as possible. Thus, t appears at most s times and the time 
evolution of GALL- is mainly determined by terms proportional to t s jt k — l/t( k ~ s ^ and hence 
GALIfc(t) - for s < k < 2N - s. 

Finally, let us consider the behavior of GALI^ when 2N — s < k < 2N. Again the fastest 
growing determinants contain the s columns of D fc with coy ^ 0, while the rest k — s columns are 
chosen among the remaining k — s time independent columns of D fc . In order to have as many 
columns proportional to t as possible these determinants should contain k—(2N—s) = k+s—2N 
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columns among the columns with i < s, as well as the corresponding r\ k columns. Thus, t 
appears at most s—(k+s — 2N) = 2N — k times and the time evolution of GALIfe is determined 
by terms proportional to t 2N ~ k /t k = l/t 2 ( k ~ N \ Summarizing, we have shown that the GALIj, 
for regular orbits lying on an s-dimensional torus behave as [22j : 

{constant if 2 < k < s 
p^r tfs<k<2N-s . (14) 
if 27V -s<k<2N 

Note that from (Tj"4"| we deduce that for s — N, GALIj. remains constant for 2 < k < N and 
decreases to zero as ~ l/t 2 ( k ~ N ^ for TV < fc < 2N in accordance with ©. 



3 Applications 

We now apply the GALI method to study chaotic, quasiperiodic and diffusive motion in multi- 
dimensional Hamiltonian systems. In particular, we consider the FPU /3-lattice of N particles 
with Hamiltonian [819) 

{q l+ i-q t ) 2 [3(q t+1 - qi) 4 ] 



with qi , . . . , qn being the displacements of the particles with respect to their equilibrium po- 
sitions, and pi, . . . ,pn the corresponding momenta. It is well known that if we define normal 
mode variables by 

i—l N 7 i—l x 7 

the unperturbed Hamiltonian (Eq. (|15[) for j3 = 0) is written as the sum of the so-called 
harmonic energies Ei having the form: 

Ei = \ {P? + ^iQf) , Ui = 2 sin ( ^W ) i = 1, • • • . N, (17) 

with LOi being the corresponding harmonic frequencies. In our study we impose fixed boundary 
conditions qo(t) = qN+i(t) = po(t) = pN+i(t) = 0, Vt and fix the number of particles to N = 8 
and the system's parameter to (3 = 1.5. 

We consider first a chaotic orbit of (|15|) . having seven positive Lyapunov exponents, which 
we compute as the limits for t — > oo of some appropriate quantities Li, i — 1, ... ,7 (see |11] for 
more details), to be a x « 0.170, cr 2 w 0.141, cr 3 w 0.114, cr 4 « 0.089, er 5 w 0.064, cr 6 « 0.042, 
0"7 w 0.020 (Fig. QJa)). We recall that chaotic orbits of Hamiltonian systems possess Lyapunov 
exponents which are real and grouped in pairs of opposite sign, with two of them being equal 
to zero. We, therefore, have in the case of the N — 8 particle FPU lattice (fT5|) a L = —an-i for 
i = 1, . . . , 8 with <7g = erg = 0. Using the above computed values as good approximations of the 
real Lyapunov exponents, we see in Figs.QJb) andQJc) that the slopes of all GALIfc indices are 
well reproduced by 1(2} . 

Turning now to the case of a regular orbit of (|15p . we plot in Fig. [2ja) the evolution of 
its harmonic energies Ei, i = 1,...,8. The harmonic energies remain practically constant, 
exhibiting some feeble oscillations, implying the regular nature of the orbit. In Figs. [2Kb) and 
[2^c) we plot the GALIs of this orbit and verify that their behavior is well approximated by the 
asymptotic formula © for TV = 8. Note that the GALI fc for fc = 2, . . . , 8 (Fig. fflh)) remain 
different from zero implying that the orbit is indeed quasiperiodic and lies on a 8-dimensional 
torus. In particular, after some initial transient time, they start oscillating around non-zero 
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Fig. 1. (a) The time evolution of quantities Li, i = 1, . . . , 7, having as limits for t — > oo the seven 
positive Lyapunov exponents cr*, i = 1, . . . , 7, for a chaotic orbit with initial conditions Qi = Q4 = 2, 
Q2 = Q 5 = 1, Qa = Qs = 0.5, Q 7 = Qs = 0.1, P» = 0, i = 1, . . . , 8 of the N = 8 particle FPU 
lattice (|15[) . The time evolution of the corresponding GALIfc is plotted in (b) for k = 2, . . . , 6 and in 
(c) for k = 7, 8, 10, 12, 14, 16. The plotted lines in (b) and (c) correspond to exponentials that follow 
the asymptotic laws © for o\ = 0.170, a 2 = 0.141, 03 = 0.114, 04 = 0.089, cr 5 = 0.064, <r 6 = 0.042, 
at = 0.020. Note that t-axis is linear and that the slope of each line is written explicitly in (b) and (c). 




Fig. 2. (a) The time evolution of harmonic energies Ei, i = 1, ... ,8, for a regular orbit with initial 
conditions qi = 52 = <?3 = qs = 0.05, q4 = gs = qe = qy = 0.1, pi = 0, i = 1, . . . , 8 of the N — 8 particle 
FPU lattice (|15[) . The time evolution of the corresponding GALI^ is plotted in (b) for k — 2, ... ,8 and 
in (c) for k = 10, 12, 14, 16. The plotted lines in (c) correspond to functions proportional to t~ 4 , t~ s , 
t~ 12 and t~ 1(i , as predicted in (|3}. 



values whose magnitude decreases with increasing k. On the other hand, the GALI^ with 
8 < k < 16 (Fig.^c)) tend to zero following power law decays in accordance to ([3]). 

From the results of Figs. Q] and [21 we conclude that the different behavior of GALIfc for 
chaotic (exponential decay) and regular orbits (non-zero values or power law decay) allows for 
a fast and clear discrimination between the two cases. Let us consider for example GALIg which 
tends exponentially to zero for chaotic orbits (Fig. [He)), while it remains small but different 
from zero in the case of regular orbits (Fig. E^b)). At f ~ 150 GALIg has values that differ 
almost 35 orders of magnitude being GALIg w 1CP 36 for the chaotic orbit, while GALIg w 1CP 1 
for the regular one. This huge difference in the values of GALIg clearly identifies the chaotic 
nature of the orbit. 
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t log(t) log(t) 



Fig. 3. (a) The time evolution of harmonic energies for a regular orbit lying on a 2-dimensional torus 
of the N = 8 particle FPU lattice (I15|l . Recurrences occur between Ei and E3, while all other harmonic 
energies remain practically zero. The time evolution of the corresponding GALI^ is plotted in (b) for 
k = 2, ... ,6,8 and in (c) for k = 9, 11, 13, 14, 16. The plotted lines in (b) and (c) correspond to precisely 
the power laws predicted in (|14[1 . 



If we select initial conditions such that only a small number of the harmonic energies Ei 
are initially excited, we observe, at small enough energies that (fT5")) exhibits the famous FPU 
recurrences, whereby energy is exchanged quasiperiodically only between the excited EiS. 

In Fig. [3l we consider the case of a regular orbit with initial conditions Qi = 2, Pi = 0, 
Qi = Pj = 0, i = 2, . . . , 8 having total energy H — 0.24. In Fig.[3fa) we see that only two normal 
modes are exited, namely E\ and £3, while all other harmonic energies remain practically zero. 
Observe that among the GALI^ plotted in Figs.^b) and (c), only GALI2 « const., indicating 
that the torus is only 2-dimensional, while all others decay by power laws whose exponents are 
the ones given in (fbl)) for N — 8 and s = 2 (see also [35] for more details). 

Let us now choose initial conditions so as to distribute the energy among 4 modes, Ei, 
i = 1,3,5,7, in our 8-particle FPU lattice. In particular, we consider an orbit with initial 
conditions = 0.1, pi = 0, i — 1,...,8, having total energy H — 0.01. What we observe 
again is that since only these 4 modes are excited (Fig. HJa)), only the GALI& for k = 2,3,4 
remain constant (Fig. [D(b)), implying that the motion lies on a 4-dimensional torus, while all 
the higher order GALIs decay by power laws, whose exponents are derived in (| 14[) for N = 8 
and s = 4. 

What happens, however, if we choose an orbit that starts near a torus but slowly drifts away 
from it, presumably through a thin chaotic layer of higher order resonances? This phenomenon 
is recognized by the GALIs, which provide early predictions that may be quite relevant for 
applications. To see this let us choose again initial conditions for our 8-particle FPU lattice, 
such that the motion appears quasiperiodic, with its energy recurring between the modes E\, 
E3, E$ and E7 (Fig. E£a)). In particular, we consider an orbit with initial conditions Q\ = 2, 
Q 7 = 0.44, Q 3 = Q4 = Q 5 = Q 6 = Q 8 = 0, Pi = 0, i = 1, . . . , 8, having total energy H = 0.71. 
This orbit, however, is not quasiperiodic, as it drifts away from the initial 4-dimensional torus, 
exciting new frequencies and sharing its energy with more modes, after about t = 20000 time 
units. This becomes evident in Fig- EKb) where we plot the evolution of Eq and E$. We see that 
these harmonic energies, which were initially zero, start having non-zero values at t sa 20000 and 
exhibit from then on small oscillations (note the different scales of ordinate axis of Figs.^a) and 
(b)), which look very regular until t « 66000. At that time the values of all harmonic energies 
change dramatically, clearly indicating the chaotic nature of the orbit. As we see in Fig. EJc), 
this type of diffusion is predicted by the exponential decay of all GALI/j , shown already at about 
t = 10000. 
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Fig. 4. (a) The time evolution of harmonic energies for a regular orbit lying on a 4-dimensional torus 
of the N = 8 particle FPU lattice (|15[) . Recurrences occur between Ei, E3, E$ and E7, while all other 
harmonic energies remain practically zero. The time evolution of the corresponding GALI^ is plotted 
in (b) for k = 2, . . . , 8 and in (c) for k = 9, 11, 13, 14, 16. The plotted lines in (b) and (c) correspond to 
the precise power laws predicted in (|14[l . 
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Fig. 5. The time evolution of harmonic energies (a) E\, £3, E5, E7 and (b) Ee, E%, for a slowly 
diffusing orbit of the N = 8 particle FPU lattice (|15[) . (c) The time evolution of the corresponding 
GALIfc for k — 2, ... ,8, clearly exhibits exponential decay already at t ~ 10000. 



4 Summary 

We have applied the Generalized Alignment Index of order k (GALIfc) as a tool for studying 
the dynamics in conservative dynamical systems, and in particular in Hamiltonian systems of 
N degrees of freedom. We have shown that these indices not only distinguish rapidly between 
chaotic and regular orbits, but also determine the dimensionality of quasiperiodic tori and detect 
slow diffusion away from quasiperiodicity, long before this becomes evident in the dynamics. 

The GALIfc represents the 'volume' of a generalized parallelogram having as edges k > 2 
initially linearly independent unit deviation vectors and are computed as the norm of the wedge 
product of these vectors. We verified numerically that for chaotic orbits GALIfc tends exponen- 
tially to zero following a rate which depends on the values of several Lyapunov exponents (see 
Eq. {2J). In the case of regular orbits lying on TV-dimensional tori, GALIfc with 2 < k < N 
eventually fluctuates around non-zero values, while for N < k < 2N, it tends to zero following 
a particular power law (see Eq. ([3|)). If, on the other hand, the orbit lies on an s-dimensional 
torus, with 1 < s < N , we showed analytically and verified numerically that GALIfc is nearly 
constant for 2 < k < s and, if s < k < 2N, decays as a power law, the precise form of which 
depends on the dimensionality of the torus (see Eq. (fT4l ). 
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Computationally, of course, it is somewhat costly to evaluate the large number of determi- 
nants needed to obtain the GALI^s by Eq. ([6|), especially in the case of large N. As a solution 
to this problem, we introduced and theoretically explained in Sect. 12.21 an efficient and accurate 
method of computing log(GALIfc) by applying the technique of Singular Value Decomposition 
(SVD) to the matrix of deviation vectors and evaluating the product of its singular values (see 
Eq. (□}). 

Applying next the GALI method to a 1-dimensional Fermi-Pasta-Ulam (FPU) lattice of 
N = 8 particles, we have demonstrated that GALIs do answer efficiently and reliably some 
fundamental questions of practical concern: When is the motion quasiperiodic and what is 
the dimensionality of the torus on which it lies? This is quite important when one wishes 
to determine the number of frequencies involved and leads to the interesting result that the 
dimension can be much less than the generally expected number of degrees of freedom N. Also, 
when is the motion not quasiperiodic but diffuses slowly away from a torus through a chaotic 
network of higher order resonances? This is notoriously difficult to ascertain (especially for 
large N), as the dynamics is 'sticky', Lyapunov exponents may be very small and thus long 
integrations are needed before the chaotic nature of the motion becomes evident. 

We believe that the results we have presented on the FPU lattice demonstrate that the 
GALIfc method can prove very useful to many practical applications: In celestial mechanics, 
there are many higher dimensional A-body problems [4123] . where one would like to locate 
invariant tori and determine the extent of phase space occupied by regular and chaotic orbits. 
In chemical dynamics, small molecules have been recently studied whose properties depend on 
the presence of stable oscillatory modes carrying large regions of quasiperiodic motion about 
them \J\. Many interesting problems are also described by iV-dimensional symplectic maps, 
as e.g. in accelerator dynamics, where the presence of large regions of tori is essential for 
maximizing the stability of betatron oscillations [5] . 

Ch. S. was supported by the Marie Curie Intra-European Fellowship No MEIF-CT-2006-025678. This 
work was partially supported by the European Social Fund (ESF), Operational Program for Educational 
and Vocational Training II (EPEAEK II) and particularly the Program PYTHAGORAS II. We would 
also like to thank both referees for their remarks. 
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